Plasmonic nanotechnology for photothermal applications – an evaluation

The application of plasmonic nanoparticles is motivated by the phenomenon of surface plasmon resonance. Owing to the tunability of optothermal properties and enhanced stability, these nanostructures show a wide range of applications in optical sensors, steam generation, water desalination, thermal energy storage, and biomedical applications such as photothermal (PT) therapy. The PT effect, that is, the conversion of absorbed light to heat by these particles, has led to thriving research regarding the utilization of plasmonic nanoparticles for a myriad of applications. The design of conventional nanomaterials for PT conversion has focussed predominantly on the manipulation of photon absorption through bandgap engineering, doping, incorporation, and modification of suitable matrix materials. Plasmonic nanomaterials offer an alternative and attractive approach in this regard, through the flexibility in the excitation of surface plasmons. Specific advantages are the considerable improved bandwidth of the absorption, a higher efficiency of photon absorption, facile tuning, as well as flexibility in the synthesis of plasmonic nanomaterials. This review of plasmonic PT (PPT) research begins with a theoretical discussion on the plasmonic properties of nanoparticles by means of the quasi-static approximation, Mie theory, Gans theory, generic simulations on common plasmonic material morphologies, and the evaluation processes of PT performance. Further, a variety of nanomaterials and material classes that have potential for PPT conversion are elucidated, such as plasmonic metals, bimetals, and metal–metal oxide nanocomposites. A detailed investigation of the essential, but often ignored, concept of thermal, chemical, and aggregation stability of nanoparticles is another part of this review. The challenges that remain, as well as prospective directions and chemistries, regarding nanomaterials for PT conversion are pondered on in the final section of the article, taking into account the specific requirements from different applications.


Introduction
With an ever-increasing demand for energy and the inevitable reduction in the dependency on fossil fuels, global energy demand looks to solar power to be a significant provider for its needs, with various solar power conversion technologies in place and rapidly progressing [1]. Electromagnetic radiation, when interacting with a material can transfer energy to its atoms, eventually converted to heat through a series of energyloss processes. This conversion of electromagnetic energy into heat is called the photothermal (PT) effect. Early stages of the PT effect were initially observed in semiconductors [2], after which researchers started to explore various material phenomena other than bandgap absorption for heat generation in nanoparticles (NPs), leading to a rapid proliferation of materials for the same. For example, organic materials undergo rapid internal relaxation by the PT effect and are often desired in cancer treatment research as they cause little damage to adjacent healthy tissues due to extremely localized heating [3]. Generally, the reduction of material dimensions to the nanoscale, such as in graphene, carbon nanotubes (CNT) and polymers, leads to an enhancement of the PT effect due to factors such as improved thermal conductivity, tunability of materials for realizing broadband energy absorption, appearance of new mechanisms of photon absorption, and improved prospects of preserving material properties [4][5][6]. Nanoparticle heating can result also due to the conversion of optical absorption by plasmons into heat. This phenomenon of surface plasmon resonance results from the interaction between electromagnetic radiation and typically high-valence materials, leading to oscillations of the free electrons in it. The decay of these collective oscillations into heat is the plasmonic photothermal (PPT) effect. The absorption characteristics such as the wavelength in plasmon resonance can be tuned and controlled by the properties of the nanoparticle such as size, shape, proximity to other particles, as well as the surrounding medium [6]. Indeed, advancements in such manipulation at the nanoscale has aided the use of plasmonic materials [7,8], such as Au nanoparticles (AuNPs), in photodynamic therapy [9][10][11]. Metal nanoparticles in general have been extensively explored in PPT applications due to their high free electron density and the possibility of intricate tuning of light absorption [12]. Noble metal nanoparticles with resonances in the UV-vis-IR part of the electromagnetic spectrum are especially researched on for PT applications [13], with excellent reviews on materials for mid-IR applications [14], cancer treatment [15], antibacterial research [16], solar-driven vapour evaporation [16], solar collectors [11,17,18], catalysis [19], clean water production [15], and wearable heaters [20,21], to name a few. This review is on PPT nanoparticle research spanning the conventional options (metals and alloys) as well as materials with induced plasmonic properties, with a special emphasis on their stability in terms of temperature and reactivity. With broad applications in therapy [22,23], laser combined imaging, solar vapour generation [24], and biosensors [25], the global market for PT devices is expected to be a multimillion dollar enterprise by 2025 [26]. This review will focus on concepts such as the theoretical aspects of PPT energy conversion (which influence material selection and design), studies on different classes and morphologies of nanomaterials that have been investigated for different applications of PT conversion, and the thermal and chemical stability of PT nanomaterials, which need to be considered prior to making the final choice. We conclude with a broad perspective on current research, challenges that remain to be solved, as well as prospects in terms of material design and deployment for better exploitation of such nanostructures for PT energy conversion.

Plasmonics in PT conversion
Of the incident radiation from the sun, 8% UV, 42.4% visible light, and 49.6% infrared radiation reach the earth's surface. Applications such as steam generation from solar power can evidently benefit from the use of materials that can absorb as much as possible of the entire spectrum of solar radiation. In this regard, plasmonic nanomaterials with tunable energy absorption can help, and the tunability of such materials only manifests at the nanoscale as changes in the absorption of incident radiation. This tunability is of utmost benefit for PT applications as the region of the electromagnetic spectrum that is not absorbed by generic PT materials can be utilized for absorption and eventual conversion into heat by incorporation of plasmonic nanoparticles of appropriate sizes and shapes. Plasmons, that is, collective electron excitations, are either excited in the bulk of the material (volume plasmons) or on the surface through excitations of the conduction electrons (surface plasmons) as shown in Figure 1. Such excitations, when occurring in nanoparticles, are termed localized surface plasmon resonances (LSPRs) as they are confined within the boundaries of the nanoparticle (in the case of continuous films, they are propagating oscillations termed surface plasmon polaritons (SPPs)). A plasmon in LSPR can be visualized as a quasiparticle confined to the volume of the nanoparticle [27]. The resulting confinement of the absorbed incident electromagnetic radiation within the nanoparticle thus means an effective localization of the incident photon energy, and the decay of this oscillation (through phenomena such as electron-electron, electron-phonon, and electron-surface scattering) releases the absorbed energy into the lattice as heat (or as photons), often making them efficient tunable PT energy materials [24,25]. Interaction of electromagnetic radiation with a material can lead to absorption, transmission, or scattering. Regarding scattering, elastic and inelastic scattering are the major classifications. Elastic scattering means conservation of the photon energy, in inelastic scattering, there are processes other than complete absorption through which photon energy can be transferred to a material. Elastic scattering is not relevant for PT applications as there is no transfer of energy into the material for heating. Absorption/inelastic scattering of electromagnetic radiation can lead to electronic, translation, vibration, and rotational transitions. The interaction time period of electromagnetic radiation with electrons is around 10 −14 to 10 −15 s. SPR falls within the  Figure 1a,b was used with permission of The Royal Society of Chemistry, from [31] ("Portable and field-deployed surface plasmon resonance and plasmonic sensors", J.-F. Masson, Analyst, vol. 145, issue 11, © Copyright 2020); permission conveyed through Copyright Clearance Center, Inc. This content is not subject to CC BY 4.0.
regime of electronic transitions and, generally, electronic transitions can be interband as well as intraband transitions. When the energy of the photons is greater than the bandgap, interband transitions are observed. As an example of the energies at which interband transitions [28] occur, Cu, Au, Ag exhibit them at 2.25, 2.4, and 4 eV, respectively, and threshold energy levels of interband transitions are 1.6-1.8 eV for Cu, Au, and Al, as well as 3.5 eV for Ag [29]. Concerning PT applications, radiative transitions such as luminescence and scattering imply inefficiency, as this scattered energy is not converted to heat. Hot electron generation and subsequent thermalization are consequences of SPR absorption that can lead to heat generation, depending on whether the decay of the SPR is through radiative or non-radiative processes. Many metals show plasmonic properties, but for PT applications there is a specific set of requirements including, but not limited to, broadband absorption of electromagnetic radiation, specifically in the UV-vis range (as infrared is already applied for heating), efficient of conversion of the absorbed energy into heat (in contrast to scattering), chemical and physical stability of the nanoparticles (e.g., against agglomeration), ease of synthesis, and low cost. Coinage metals, such as Au, Ag, and Cu, with high densities of free electrons exhibit plasmon resonances in the visible region suitable for PT applications [30]. The subsequent part of this review highlights the parameters that influence various properties of plasmonic materials relevant to PT energy conversion. We first derive the equation quantifying the absorption frequency of plasmons, followed by a discussion on the changes to this frequency that can be induced by changing the governing parameters of this equation, and conclude with a few examples that model the optical scattering properties of generic morphologies such as spheres and nanorods.

The plasmonic oscillation frequency
The optical response of plasmonic nanoparticles, such as AuNPs, to incident electromagnetic radiation depends on their size, shape, morphology, proximity to one another, as well as the surrounding medium [32]. The vast changes in absorbance due to changes morphology stem predominantly from changes to the directionality of the LSPR (due to changing curvature and dimensionality of the nanoparticle) and changes to the dynamics between the restoring and exciting forces of the plasmons, such as the mean free path, the relative contributions to the plasmon damping of different scattering phenomena, the different scattering processes of the oscillating plasmons, and the screening between the plasmons and the restoring nuclear forces [33,34]. The proximity among LSPR-active nanoparticles is also a major factor. Indeed, combined effects of proximity as well as morphology influence considerably the LSPR properties, for example, in Au nanorods and nanospheres. In contrast, nanospheres and nanorods exhibited considerable tunability of the LSPR due to changes to the localized electromagnetic field of the plasmons due to changing curvature [35]. Finally, changes to the material composition, such as through doping or vacancy processing, can affect the LSPR because of changes in the free electron density, the electron effective mass, and the electronic band structure in general [36,37]. An understanding of the changes in absorbance with respect to changing parameters of the material under consideration can be developed with a few examples of the theoretically arrived optical cross sections of a few generic morphologies of nanoparticles and will be discussed next for the case of nanospheres, nanorods, and nanomatryushkas.
The arrival at the expression for the LSPR frequency of a free electron cloud (as is typically assumed to be present in metals) starts with the relations between the dielectric displacement (D) of the electron gas in relation to the incident electric field (E) which it is [38] subjected to, given by (1) wherein P is the polarization density. P can be arrived at by solving the equation of motion for a single electron as (2) Hence the expression relating the dielectric displacement (D) and the external electric field can be obtained as (3) where is the natural frequency of oscillation of the electron cloud. Comparing Equation 3 to the general constitutive relation for a linear isotropic material given by Equation 4, we get the relation in Equation 5. ε r is the relative permittivity of the material and ε 0 the permittivity of free space.
For frequencies close to ω p , the temporal duration of damping (quantified by the product of ωτ, where τ is the relaxation time of the free electron gas) is much higher than unity, thus leading to an approximation that there is no damping. Hence ignoring the damping term in Equation 5, we get (6) It follows also that under plasmon resonance conditions ε 1 < −ε m , ε m is the dielectric constant of the surrounding medium and hence the LSPR frequency can be arrived at as [39]: Thus, for any morphology of a plasmonic nanoparticle, the LSPR frequency is intimately tied to the free electron density and the dielectric constant of the surrounding matrix. These factors thus decide the shape, position, and width of the plasmonic absorption and will be further elaborated in the subsequent sections.
Plasmon absorption is also determined by the nanoparticle shape, which, although it does not appear in the plasmonic frequency equation, manifests as a shape/size factor in the calculations of the extinction spectra of nanostructures. Illustrative examples of the absorbance spectra for different morphologies of Ag nanoparticles are shown in Figure 2, elucidating the influence of the same. The extinction spectra (the summation of absorption and scattering spectra) from which the shape effects of different morphologies on plasmon excitation can be understood are hence crucial for assessing the PT properties of nanomaterials.

Tunability of the plasmon frequency -changes to the dielectric constant.
In 1834, William Whewell coined the term dielectric [41]. In a dielectric material, positive charges are arranged in the direction of electric field and negative charges opposite to the field, causing a polarisation. The wavelength-dependent electric field and dipole moment determine the dielectric property of the material and also the plasmon absorption (by affecting the polarizability, as was shown in prior sections) [42].
The real (ε r ) and imaginary part (ε i ) of the dielectric constant relate to the refractive index (n) and extinction coefficient (k), respectively. Energy from a time-varying incident electric field is dissipated in part as heat, termed as dielectric loss. It can be envisioned that this dielectric loss is hence an important attribute to be considered for PT applications. The dielectric loss (δ) can be expressed as: (8) where ε i is the imaginary part of the dielectric function, and ε r is the real part of the dielectric function [43]. Polarization occurs for charge carriers in a number of orbitals or band, which then very often overlap, leading to multiple transitions even during plasmonic resonance. For example, metals such as Au, Ag, and Cu have d-band electrons close to the Fermi surface, and the polarization (P) in the presence of an electric field (i.e., the plasmon oscillation) has contributions from interband transitions as well. The dielectric function must account for this, and thereby the conventional definition given as [44]:  changes (for metals at optical frequencies, for example) to (10) χ ∞ is the susceptibility arising from the core electron polarizability (causing interband transitions), and χ D is the corresponding susceptibility of the conduction electrons (modelled through the Drude assumption of a free electron "sea"). Hence, PT applications targeting plasmonic materials must account for the contributions to the dielectric function of the terms mentioned in Equation 10. A broad absorption of wavelengths is a reinforcing attribute of plasmonic materials for PT applications. However, increased broadening is associated with a reduced absorbance due to an increased scattering of the plasmon oscillations. An important phenomenon that broadens the plasmon line width is the scattering of plasmons at the surface for nanoparticles of sizes approaching the mean free path of electrons. Apart from the minor contribution to the linewidth arising from the disparity in particle sizes (when considering the absorption of a cluster of nanoparticles), the linewidth is controlled solely by surface scattering in such nanoparticles [45]. For free electron metals, the frequency of electron scattering (γ) is equal to the width of the plasmon resonance (Γ) [46]: (11) In order to account for surface scattering, the linewidth needs to be modified as, (12) where ν F is the Fermi velocity of electrons, Γ 0 is the plasmon linewidth of the particles or damping constant, a is the radius of the spherical metal particles, and A is a parameter that depends on the scattering process. The final expression for the dielectric function taking into account the inter-and intraband transitions as well as the changed linewidth due to surface scattering (as appropriate), is (13) Multiple inferences can be made from this formulation. In addition to linewidth broadening contributions to the dielectric func- Figure 3: Manifestation of the governing factors of SPR. Shown are the changes to the peak position while considering only the shape (red circles), the shape as well as the electron density (blue squares), and the combination of shape, electron density, and deformation potential (green diamonds). Figure 3 was reprinted with permission from [53], Copyright 2017 American Chemical Society. This content is not subject to CC BY 4.0. tion for smaller nanoparticles, electron orbit contractions (due to the majority of electrons being in proximity to the surface) result in an increased Coulombic force of restoration and hence a shift in the dielectric function [47]. Similarly, with a decrease in grain size, due to the fact that there is an increase in the volume fraction of grain boundaries compared to the grains, and since the dielectric strength of a grain is lower than a grain boundary, the dielectric permittivity decreases with decreasing grain size [48]. Moreover, the interaction between plasmonic nanoparticles and substrates on which they are deposited cannot be ignored. The polarization of charges in the nanoparticles induces dipoles in the substrate atoms in proximity of this polarization field, which in turn affects the nanoparticle resonance. This has been observed to induce higher-order resonances when the mismatch between the permittivity of the substrate and the surrounding of the nanoparticle increases, as well as to (depending on the orientation of the applied field with respect to the induced fields in the nanoparticle and the substrate) increase or decrease in the absorption intensity [49][50][51].
It follows from the discussion that considerable effects to the plasmon resonance are affected by the dielectric surrounding the plasmonic nanoparticle, apart from the permittivity of the nanoparticle itself. The shift in resonance on incorporation of the nanoparticle into a dielectric will decide its absorbance and hence the efficiency of conversion of the resonance into heat. This has been verified experimentally in multiple studies, wherein plasmon absorption peak shifts of up to 150 nm [52] have been observed with just a unit change in refractive index of the surrounding medium. As important as the dielectric function, the free electron density of the nanoparticle is also crucial for PT applications. This factor decides multiple attributes of a material, such as the quantum of this value that makes a given material plasmonic, the shifts in absorption and linewidth, and the ratio of radiative to non-radiative to radiative damping. A discussion on the manifestation of the free electron density on the plasmonic performance of different nanomaterials follows and is illustrated in Figure 3. 2.1.2 Tunability of the plasmon oscillation frequencychanges to the free carrier density: Free carrier densities of electrons are in the range of 10 22 to 10 23 cm −3 for plasmonic metals such as Au, Ag and Cu [54]. The free electron density is tied to the effective mass and determines the resonant plasma frequency (ω p ), given by (14) where N is the free electron density, m is the effective mass, and e is the elementary charge. The LSPR effect is not present in most of the semiconductors because of their lack of the required free carrier concentration. Similar to how the free carrier density of metals can be tuned by size, morphology, and refractive index of the nanomaterial, the free carrier density of semiconductors can be easily tuned by doping, temperature variations, or by phase transitions. LSPR in semiconductors in the NIR-mid-IR region is possible when the free carrier concentrations lies between 10 16 and 10 19 cm −3 [54].
As an example of how the free electron density influences the plasmon resonance when materials with different work functions are combined, metals in contact with semiconductor metal oxide nanoparticles exhibit the spill-over effect, which alters the plasmonic absorption and spectral width of the plasmonic nanoparticles integrated in dielectric matrices. This spill-over effect, however, decreases with a decrease in electron density. For nanoparticles with low electron density (typically for radii less than 10 nm), the omnidirectional diffuse scattering dominates the resonance and the spill-over effect can be safely neglected. The spill-over effect and diffusive scattering can be related by [55]: (15) m* is the effective mass of an electron, r s is the radius of on-electron, d T is the complex length, is the classical Mie LSPR frequency of the sphere = , and the Bohr radius is The free electron density can be tuned to make non-plasmonic materials plasmonic, which is useful for PT applications that require materials with properties that are not present in conventional plasmonic materials, such as a higher melting point, alloying capabilities, as well as possibly lower reactivities. Morphologies with optimal PT properties in terms of optical absorbance can benefit well from prior knowledge of the same, which can be obtained quite well through modelling efforts. Modelling of the absorbance of a plasmonic nanoparticle is done by calculation of its optical cross sections and specifically the extinction cross section, which includes the absorption as well as the scattering cross section. A few examples of the same will be discussed to compare the differences that need to be accounted for regarding different morphologies, essential for obtaining the accurate results in Figure 4.  Figure 4 was reprinted with permission from [56], Copyright 2015 by the American Physical Society. This content is not subject to CC BY 4.0.

Nanospheres:
The interaction of light with a particle is in one way simplistically modelled using the so-called quasi-static approximation where the incident electric field is assumed to be spatially uniform. This assumption is valid only for wavelengths much larger than the particle size. The extinction cross section, which is the result of this modelling and which is an expression summing up absorption and scattering of the incident radiation, is derived starting from Laplace's equation with an electric potential (ϕ), (16) Calculation of the resulting scalar potentials inside and outside the particle leads to the expression for the polarizability (α) of the particle, (17) ε m is the dielectric constant of the surrounding medium. The extinction coefficient (σ ext ) is introduced as: where C ext is the extinction cross section. On substituting for the polarizability, the final expression for the extinction cross section [57] for a spherical particle interacting with light is arrived at as (20) where a and ε m are the particle size and dielectric constant of the surrounding medium, respectively, and ε 1 and ε 2 are the real and imaginary parts of the dielectric function of the material.
For larger particles (above approx. 40 nm), field-retardation effects affect the resonance position due to increasingly higher radiative damping, with significantly different extinction spectra arising also from the excitation of multipole resonances, which is not captured by the quasi-static approximation as shown in Figure 5.
The Mie theory solution of the extinction cross section is used to account for this through the inclusion of a size parameter (x = 2πa/λ) as well as through a coefficient for covering partial Figure 5: Extinction efficiencies of gold nanospheres calculated through the quasi-static approximation vs using Mie theory. The expected trends in the plasmonic peak shifts with particle size post the quasi-static limit can be seen to be predicted better by Mie theory. electric and magnetic waves (for the multipole orders). The extinction cross section is then given by (21) where Here, ψ and ξ j are the Riccatti-Bessel functions, and J and Y are Bessel functions of the first and second order, respectively, m is the ratio of n m and n, the real refractive index of the surrounding medium and the complex refractive index of the spherical nanoparticles, and j is an integer representing the order of scattering (dipole, quadrupole and so on). These equations for the extinction cross section more accurately capture the absorbance for particles larger than 40 nm. Even by accounting only for the dipole mode in such particles, depending on the order to which ψ j and ξ j are expanded, the shift as well as the changes to the peak broadening can be captured. For example, expanding the function to the order of x 2 yields [58]: (22) This equation captures well the redshift in the plasmon resonance for increasing particle diameters, not observable with the quasi-static approximation alone, as shown in Figure 6 for the case of Au nanospheres. Redshift and broadening of the plasmon resonance for the simple spherical morphology, allowing tunability of the absorption characteristics, are also evident from this discussion.
A common analogue to the spherical morphology are concentric nanoparticles, wherein a core nanoparticle is surrounded by a shell layer of specified thickness and composition. The change in refractive index of the shell in comparison to the core leads to a multitude of interesting effects such as altered/enhanced absorption and stability [60]. Hence core-shell nanoparticles can be tuned very effectively to the desired wavelength range by manipulating the thickness and/or composition of the shell in addition to the tunability prospects of the core itself, such as size and composition.

Nanorods:
Nanorods in the shapes of cylinders or rectangles are a very common morphology explored in applications of sensing, catalysis, and plasmonics. Their attributes, including aspect ratio (ratio of the length to width), curvature, homogeneity, dimer formation, and placement, make them very interesting for multiple applications such as communication, hotcarrier enhanced catalysis, high-temperature sensing [61], and electronics applications such as transistors. Their extinction coefficient is often calculated using Gans theory for randomly oriented nanorods (with the rod geometry assumed as prolate spheroids with three principal axes) in the dipole approximation: (23) where ε m is the dielectric constant of the surrounding medium, P j is the depolarization factor along the A, B, and C axes, and R is the aspect ratio of the nanorods (B/A or width/length). For example: Because of the distinct and significant changes to the polarization of the electron cloud with respect to the elongation (aspect ratio) of the rods and any such rod-like nanostructures, considerable changes to the absorption properties are observed. This is because in addition to the polarization changes, there are also changes in the scattering processes (radiative vs non-radiative) of the plasmon oscillations. It follows that elongated nanostructures of plasmonic nanoparticles are more conducive for PT applications than spherical particles due to the precise and disparate manipulations possible in the absorbance of the former.

Matryushka:
Multilayered nanoparticles have interestingly different optical characteristics than their single or nonlayered pristine counterparts. Taking multiple forms such as core-shell, sandwich structures, and films, such multilayered nanostructures exhibit tunable plasmon resonance bands, resulting in additional shifts due to plasmon hybridization [62]. As an example, the interaction and hybridization of the plasmons of two separate metal shells (nanomatryushka) and two core dielectrics constitute the nanomatryushka (simplistically viewed as a ring within a ring). The coupling strength and energy between the plasmons on the inner and outer shells will determine the plasmon properties, which can be tuned by altering the dielectric spacer layer as well as the thickness of the specific nanoshells. Because of this intensified interaction between inner-core and outer-shell plasmons, PT energy transduction is significantly more effective.
The absorption coefficient (C abs ) of the nanomatryushka (NM) can be calculated by varying the volumetric factor of the different layers of the nanostructures and the refractive index of the surrounding medium [63]: (24) P LSPR is the power absorbed by the nanoparticle, ε 0 and μ 0 are, respectively, the permittivity and permeability of vacuum, n m is the refractive index of the surrounding medium, and is the amplitude of the transverse incident wave of the laser source. Equation 24 can be simplified in the quasi-static limit as: (25) with h is the height of the nanoshells, R a , R b , R c , and R d are the radii of the different concentric spheres, λ LSPR is the LSPR wavelength of the nanoparticle, and ε is the dielectric constant of the surrounding medium. In this nanomatryushka morphology, h represents the distance between the two concentric spheres, whereas in spherical nanoparticles the parameter inducing size effects was different. Hence, it is conclusive that the quasi-static approximation is valid for different morphologies but needs to account for the size when exceeding the validity of this approximation. Simulated absorption cross sections of Ag nanomatryushkas for different radii of different concentric spheres are shown in Figure 7.  From the few examples presented, it is clear that information such as the absorption range, the intensity of absorption, and importantly the changes in absorption resulting from changes to shape, size, material, and dielectric can be effectively obtained from simulations. Such endeavours minimize the effort required for material optimization. Also, they allow for considerable flexibility in exploring unconventional morphologies and material combinations. With continually improving computational resources, allowing for more complex calculations of extinction properties as well as PT conversion, modelling techniques are expected to play a symbiotic part with experiments in PT research. Understanding the processes that lead to the conversion of light to heat is critical for predicting and understanding PT performance, for which a discussion on the same is presented next.

Relaxation mechanisms of plasmons -the conversion of photons to heat
When photons strike a plasmonic nanoparticle, plasmon excitations result in changes to the electron cloud thermal distribution. The electromagnetic absorptivity of the metal conduction electrons is proportional to the rate of momentum transfer to the lattice. It can be written as [64]: (26) m * is the effective electron mass, n is the electron density, and τ eff is the conduction electron relaxation time given by (27) τ is the conductivity relaxation time with temperatures T larger than the Debye temperature Θ. The Debye temperature is the temperature of a crystal's highest mode of vibration. The decay of the excited electrons (plasmons) is through either radiative relaxation (i.e., photon emission) or non-radiative relaxation. Non-radiative relaxation occurs through the sequential processes of hot carrier generation, electron thermalization, and finally electron-phonon coupling heat transfer [65], leading to a temperature increase of the lattice ( Figure 8 [66] and Figure 9). Thus, the energy exchange between energetic plasmons and the phonon modes leads to an increase in temperature of the metal nanoparticles, and the ultimate conduction of this heat to the surrounding. The exchange time will vary for small nanoparticles due to quantum confinement effects [65] observed when the size of the nanoparticle is comparable to the Bohr exciton radius. In such cases, the confined electron wave functions have discrete energies and undergo considerably higher scattering with the surface because of the "shrunken" wave functions, contributing significantly to the reduction in absorption due to significantly increased linewidth broadening [68]. Interestingly, such small nanoparticles (typically smaller than 5 nm) have been found to undergo a much higher increase in temperature than larger nanoparticles (by more than five orders of magnitude) on irradiation with a single photon due to a much more efficient energy conversion [69,70]. Inferences from studies on PT relaxation in plasmonic nanomaterials are useful in designing them for such applications.  Electron-electron scattering can be understood to be the process following hot carrier generation due to plasmon relaxation, wherein thermalization of the electron distribution occurs. For a single scattering event, the duration of electron-electron thermalization is generally of the order of 10 fs (depending on the starting energy E) [52]. Accordingly, the influence of the initial energy on the electron-electron scattering time (τ e-e ) with initial energy (E) and a Fermi energy (E F ) can be expressed as (28) Electron-electron scattering is followed by electron-phonon coupling and heat transfer to the surrounding [71,72]. These non-radiative decay processes specifically concern PT applications.
Electron-phonon coupling, which reduces with decreasing particle size, and electron-surface scattering, which rises with decreasing particle size define the relaxation dynamics of the photoexcited electron gas into heat [52]. Interestingly, for Ag and Au, the transduction of absorbed light into heat due to electron-phonon coupling was found to be independent of the particle size above 10 nm diameter (though surface scattering effects start to dominate when the particle size is reduced beyond the mean free path of 40 nm for these nanoparticles) [73]. Electron-electron thermalization timescales for noble metals such as Au range from 10 to 100 fs, whereas the time scales for the electron-phonon interaction are slightly shorter, such as 10-30 fs for Al [74]. Hot carrier lifetime and electron-phonon coupling of selected metals are shown in Figure 10 and Figure 11. From the representation of electron filling for the different materials shown in Figure 10, it can be seen that the carrier lifetimes (depicted by the curvature of the  Figure 11a and Figure 11b were reprinted with permission from [75], Copyright 2003 by the American Physical Society. This content is not subject to CC BY 4.0. Figure 11, right panel was reprinted with permission from [76] , Copyright 2000 by the American Physical Society. This content is not subject to CC BY 4.0. Fermi surface) can vary considerably due to differences in the scattering rates of electrons that occupy the different energy levels.
The differences in the mean free path of carriers with the difference in their energy from the Fermi energy validates the observation that the higher the difference of the electron energy from the Fermi energy, the lower the mean free path. Therefore, it is essential for an understanding of the distribution of electrons based on their energy to predict the carrier lifetimes of plasmons, which in turn decides the damping properties of the resonance.
Since all mechanisms of plasmon relaxation (through non-radiative scattering) result in the eventual coupling of the scattered plasmons to phonons, their understanding is critical for designing materials for such applications. Simplistically, phonons are quantized particles or waves with a unit of vibrational energy corresponding to the continuum of frequencies ω(q) [77]. Their postulation and confirmation arose initially to reconcile deviations of the specific heat of materials at different temperatures from the Dulong-Petit law. The quantized characteristic (Eigen)energies of these vibrations are simply:   where ℏ is Planck's constant, and q is the wave vector of the phonon oscillations. These elementary excitation phonons are bosons with a wave vector of k = 2π/λ. Crystal defects and other such sources of anharmonicity can change their frequencies and, hence, their coupling characteristics to electrons and the ultimate conversion of electron scattering into heat.
The energy distribution of hot carriers (which decides the relaxation times) depends on the electronic band structure [78], particle size, density of states, and the geometry of nanoparticles [79]. Figure 12 shows the vast differences in the population of electrons and holes for Al, Ag, Au, and Cu. This distribution takes the characteristic of a symmetric/continuous distribution of electrons and holes in aluminium, a bimodal distribution in silver, and an antisymmetric hole-dominant distribution in copper and gold.
Also, decreasing the thickness of Au films from 40 nm increases hot carrier generation through increased intraband transitions as well as due to an asymmetric density of states distribution. Hence, hot carrier generation and lifetimes depend on the material, particle size and shape, as well as the distribution of the density of states [81]. Figure 13 illustrates the relaxation process of the phonon vibration for a better understanding of phonon dynamics.
The electron-phonon coupling constant, which describes the potential of a material for undergoing lattice heating, is calculated by transient reflectivity studies, for example through pump-probe laser spectroscopy [83]. From the slope of this transient reflectivity curve, the coupling constant is arrived at, from which the lattice temperature is calculated through the two-temperature model by means of two coupled diffusion equations, where T e and T l are, respectively, the electronic and lattice temperature, C e and C l are, respectively, the electronic and lattice heat capacities, g is the electron-phonon coupling constant, and T 0 is the ambient temperature.
It can be seen that the governing fundamental particles in PPT heating are phonons, and that the tailoring of PT heating properties by phonons can thus be effectively realized by lattice tuning of the pristine materials, that is, by defect tailoring, doping, or by providing external stimuli such as pressure.
While considering the actual increases in temperature of the nanoparticle possible due to plasmonic heating, the differences in specific heat capacity between the lattice and an electron must be considered. Electrons having a much smaller specific heat capacity compared to the lattice will undergo much higher heating than the lattice as a whole. For example, for differences in the electron temperature (due to the initial thermalization step) up to 1000 °C, the increase in nanoparticle temperature was found to be only 20 °C [84]. This implies by no means a weak efficiency of heat production due to plasmon relaxation, as high temperatures (230 °C for AgNPs for example) are possible, depending on the various factors that have been discussed in preceding sections [85]. PPT heating has subwavelength spatial selectivity, implying that changes to absorption and, hence, heating can be affected by means of changing the lasing power density and/or the wavelength. Indeed, NP morphologies have even been shown to be controlled due to extremely localised heating [38].
For increasing the spatial distribution of the generated heat, compared to an individual nanoparticle, nanoparticle assemblies and/or increased interfaces (by including constructions such as holes and other scattering centres) within a single nanostructure can be effective [86][87][88][89]. A few governing equations which allow to predict the different aspects of PPT energy conversion are discussed next. They can be of use to explicitly select a material/morphology/matrix based on the application it is being considered for.

Calculations in PT heating
Calculations of PT heating metrics, such as the specific heat capacity of materials with added nanoparticles, the conversion efficiency, the rate of heat generated per volume, and the temperatures at different locations across the thermal front require appropriate equations. The PT conversion efficiency for water vaporisation can be calculated from (32) Δm is the water mass loss during irradiation, is the phase change enthalpy of water to vapour, M is the molar mass of water, I is the solar power intensity at the surface of the sample, S is the irradiated area of the water surface, and τ is the irradiation time. This equation is applicable for the case of nanoparticles dispersed in liquids for vapour generation, the most common non-medical application of PT energy conversion [90]. The equation governing the temperature distribution in an isotropic nanomaterial is given by (33) where ρ is the density, C is the specific heat, k is the thermal conductivity, T is the temperature, t is the time, and is the rate of heat generation (or depletion) per unit volume. The Laplacian ∇ 2 represents the temporal change in the temperature variation of the material. The rate of heat generation or depletion, , is arrived at from [91]: (34) Solving the second-order Equation 33 requires specification of boundary conditions as well as initial conditions such as the type of heat transfer applicable to the specific study (conduction/convection/radiation), the geometry (radii, lengths, and similar dimensionally defined factors), direction of heat transfer, and the heat transfer coefficients for the material as well as the surrounding medium [84]. Thermal capacitance coefficient (β) for the ellipsoidal, rod, disk, and ring morphologies of Au nanoparticles as a function of the aspect ratio D/d. Figure 14 was reprinted with permission from [96], Copyright 2010 American Chemical Society. This content is not subject to CC BY 4.0.
As the method of laser irradiance is generally employed in PT studies, to measure the temperature rise of the nanoparticles when irradiated with a laser source, the energy balance equation can be used [92,93]: (35) m i and C i are, respectively, the mass and the specific heat capacity of the i-th component, T and t are, respectively, temperature and time, and Q in and Q out are, respectively, the rate of input power from the laser and the rate of heat dissipation to the surrounding. Often, nanomaterials are used in conjunction with a phase-change material for energy storage applications, and when plasmonic nanoparticles are integrated into a solid phasechange material (n-PCM), the energy balance equation is be given by [94]: (36) ΔH fus is the heat of fusion, ρ s is the density of solid, ϕ is the volume fraction of the particles, T m is the melting temperature, T c is the temperature of the cooling fluid, W is the thickness of the compartment holding the n-PCM, k w is the conductivity of the wall, and k′ is the conductivity of the PCM. The thermal conductivity of the plasmonically enhanced PCM materials can be varied with the volumetric fraction of the PCM material and can be specified by [95]: , K np and K PCM are the thermal conductivities of nanoparticles and PCM, respectively, ϕ is the volumetric fraction of the nanoparticles, ρ np and ρ PCM are the densities of nanoparticle and PCM, respectively, and C pPCM is the specific heat capacity of the PCM.
The thermal capacitance coefficient relates the rise in temperature of the irradiated nanoparticles to the thermal conductivity and absorption power of the surrounding medium. The temperature increase in the metal nanoparticles (ΔT NP ) and thermal capacitance coefficient can be calculated using the following equation [96]: (38) I is the illumination intensity of the laser source, σ abs is the absorption cross section of the nanoparticles for the relevant wavelength, β is the thermal capacitance coefficient, R is the radius of the sphere, and k med the effective permittivity of the homogeneous medium. The thermal capacitance coefficients (C th ) of the different morphologies mentioned in Figure 14 can be calculated as: (39) Employing the appropriate of the described equations ensures compatibility of assumptions and the observed results. Any deviations could imply a multitude of factors, such as differences in assumptions in terms of the dominant mechanism of heat transfer, a deficient representation of the particle geometry, or inaccuracies in modelling the temporal evolution of the generated heat, to name a few. With this understanding of the myriad of theoretical and phenomenological phenomena, considerations, and evaluations in PPT heating, the different materials that have been investigated are well worth reviewing to appre-ciate the intricacies in this field. This will be done while keeping in mind two important factors. One will be the different materials that are involved in PT heating, and the other will be the stability requirements for the materials used and possible means to improve this stability.
4 Plasmonically active nanomaterials for PT applications 4.1 PT performance of metal nanoparticles Metallic nanoparticles are a dominant group of elements for PT research. The calculations of different optical cross sections, namely absorption, scattering, and extinction allow for the prediction of the ideal candidate for maximizing PT efficiency. For example, in terms of maximizing the absorption and, hence, the interaction with solar radiation of these nanoparticles for steam generation, modelling of these cross sections with water as the surrounding medium can give precise results [97]. The ratio between absorption and scattering cross sections can be intuitively concluded to be a convenient measure of the absorption efficiency, as maximizing absorption whilst minimizing scattering implies a higher conversion of light into heat rather than re-emission as light. Considering the broad spectrum of metallic elements, nanospheres of gold, silver, platinum, cobalt, zinc, nickel, titanium, copper, aluminium, molybdenum, vanadium, and palladium have been analysed [98]. Although there were assumptions of homogeneity, sphericity, and no interaction between particles, the Mie theory-based scattering cross sections in this study can still be important for the assessment of optical properties. In this study the conclusion that all particles become increasing scatterers with size (from 10 to 50 nm), with an analogous trend in the maximal temperature increase of the nanoparticles, was made. The larger nanospheres (larger than 50 nm) acted more as absorbers than as scatterers for wavelengths of 200-500 nm for Au, 400-1000 nm for Pt, 600-1000 nm for Pd, Co, and Zn, 500-1000 nm for Mo, 200-600 nm for Cu, 300-1000 nm for Ni, 450-1000 nm for V, and 200-1000 nm for Ti. Al and Ag were inefficient absorbers and predominantly scattered the light for radii between 200 and 1200 nm. A detailed investigation on Au, Ag, Cu, and Al with regard to solar heating has been done wherein the effects of changing size and morphology of the nanoparticles and the dielectric environment were studied. An intuitive but important factor for selecting PT materials is the ratio between the absorption and scattering sections, termed the absorption ratio (α) and given by: I(λ) is the intensity of the incident radiation, C abs (λ) is the absorption cross section, C scat (λ) is the scattering cross section, and λ max and λ min are the maximum and minimum wavelength of the calculated solar radiation, respectively. The higher the ratio, the more suitable the material is for PT heating. With this in mind, a few of the key insights are the highly intense but narrow absorption cross section [99] of Au compared to other metals, the higher δ of Al and Cu, the higher temperature increases possible with reduced sphericity, the better absorption properties when surrounded by materials with a higher refractive index (due to multipolar excitations instead of only dipolar excitations), and the threshold particle size above which the optical scattering starts to increasingly manifest, reducing the PT heating efficiency (e.g., 30 nm spheres of Au). Such calculations can help in choosing the right plasmonic metal for an application, and since PT applications rely on maximal absorption by the plasmonic resonance, information on the required wavelength range and the available particle sizes can narrow down the choices.
Metal nanoparticles in general are very efficiency in generating plasmons since they have many mobile electrons with a mean free path of 10 to 100 nm [100]. They are exceptional in their tunability of wavelength absorption due to tunability of shape, from 0-dimensional to 3-dimensional nanoparticles, size, and composition. The coinage metals Ag and Au have been explored traditionally for PT applications. A schematic representation of the plasmonic absorbance regimes of Au, Ag, and Cu for radii between 50 and 100 nm is shown in Figure 15. Au nanoparticles have been extensively used for PT applications because of the possible intricate control of morphology (and hence absorbance), the chemical stability, and the obvious plasmon-enhanced heat transfer that does not occur in low-cost, yet ineffective, non-plasmonic materials such as graphite and CNTs [101][102][103].
Excellent reviews on properties relevant to PT heating such as the relaxation dynamics, field distributions during SPR, shape and volume effects, and the optimal configuration of Au and Au-based nanomaterials in PT heating systems further support this view [87,104]. Careful studies done on plasmonic Au nanoparticles revealed the different mechanisms that can lead to vapour formation in steam generation. Particularly, the influence of the irradiance intensity on two contrasting mechanisms, namely explosive boiling and plasmonic nanobubble formation has been revealed [83]. Plasmonic nanobubbles (formed when the irradiation fluence exceeded a threshold value), although being excellent tunable scatterers themselves, did not result in thermal phenomena such as heating and only led to mechanical phenomena such as cavitation effects. Explosive boiling is of explicit use in PT applications, while its use in a novel synthesis strategy for SERS-enhancing self-assembled nanoparticles is worth mentioning [105]. Though details about the absence of thermal phenomena for plasmonic nanobubble formation are yet to be confirmed, the important conclusion that increasing the irradiation fluence will not necessarily lead to better PT properties can be made from such studies. This can be understood from a very helpful phase map ( Figure 16) [106] showing the different phenomena that accompany the laser irradiation of plasmonic nanoparticles. The influence of irradiation fluence, pulse duration, continuous vs pulsed irradiation, morphology, and solvent has been studied in detail for Au nanoparticles regarding applications such as catalysis, bio-activation, steam generation, and thermoplasmonics.
For PT applications, polygonal Au structures are preferred over nanospheres because absorption wavelengths can be tailored and broadband absorbance is possible, in contrast compared to a single narrow absorbance peak for nanospheres as discussed in preceding sections. A ubiquitous example of such polygonal morphologies, the Au nanorod, has been investigated quite well with regard to the effect of the aspect ratio as well as different encapsulating shell materials (Ag 2 S and ZnS) on PT perfor-mance. The temperature increase depends on the difference between the exciting wavelength and the plasmonic wavelength, with minimal differences contributing to maximal temperature increase. An associated increase in the scattering cross section calculated theoretically validated this observation as scattering increases when the difference between the exciting wavelength and the plasmon wavelength increases. This difference can be reduced by shifting the resonance wavelength closer to the excitation wavelength, possible through incorporating a shell that can cause this shift due to changes in the medium refractive index, as can be understood from Equation 22. Shell materials can also increase PT efficiency by absorbing radiation over a broader range inaccessible by the plasmonic core [89]. The separation between excitation and resonance wavelength can also be reduced by inducing the assembly of isolated nanoparticles, such as the assembly of Au rods and polyhedral particles in the presence of glutathione.
Reports on the relative radiative and non-radiative damping contributions to plasmon decay in nanorods and nanospheres revealed that the magnitude of the scattering cross section (and hence PT conversion) depends on size and shape rather than the volume, with smaller particles being efficient absorbers and larger nanoparticles effective scatterers [92]. The comparison between Ag, Pt, and Pd confirmed this. Non-radiative decay was absent in Ag above a certain particle size, while Pt and Pd exhibited non-radiative decay for all sizes, confirmed through conductivity measurements on plasmonically generated hot electrons. Another comparison between plasmonic heat generation in Au nanorods and spheres revealed a much higher (60%) heating efficiency for nanorods, due to the presence of more atoms near the surface (due to the flattened morphology of rods compared to spheres, Figure 17) which contribute more to the heating effect than the interior atoms that do not participate in the resonance [107]. Such improvements to PT heating because of a higher surface-to-volume ratio has also been studied with shapes such as rings, ellipsoids, rods, and disks [96], indicating another possible approach to maximize PT efficiency.
Remarkable changes to absorbance manifest also through changing the proximity between nanoparticles. Since the concentrated electric field distributions during SPR extend well outside the physical boundaries of the plasmonic structures, their influence on adjacent particles can lead to interesting resonance phenomena. As an example of such interactions, the most common arrangement of a nanoparticle dimer is often considered. In general, the interactions are classified into repulsive and attractive, wherein the former blueshifts the resonance with an associated increase in damping. In contrast, attractive repulsions lead to an increase in resonance intensity coupled with a redshift. Multiple phenomena come into play while considering the coupling between nanoparticles, such as separation distances, morphologies, and the influences of dipole-dipole/dipole-multipole/multipole-multipole coupling, to name a few. Generating a universal model that describes plasmon coupling is not possible. As an example of this, the conventionally used dipole approximation, which predicts the decay of the near field fails for small interparticle distances (at separation-to-diameter ratios less than 0.1). For such cases, alternate models explain the observed results, such as the plasmon hybridization model and extended multipole methods. The plasmon hybridization model, for example, considers the interactions at such distances through the formation of bonding and anti-bonding plasmon modes, yielding a robust prediction of the resonance shifts [108,109]. As the particles approach, more modes occur, indicating that the system's symmetry reduces to a side-to-side approximation. When two nanoparticles interact, their unique modes pair to form new dipole resonances of different energies. For example, in nanorods arranged end-to-end, the net dipole moment of the higher energy mode is zero. In an asymmetric nanorod dimer (heterodimer), both bonding and anti-bonding plasmon modes occur at the same time, resulting in a non-zero dipole. The fractional shift in resonance wavelength depends on the size (S), the separation between the nanoparticles (R), and shape factors (Λ, γ): A comprehensive treatise of the fascinating topic of plasmon coupling is outside the scope of this review, and hence only a rough overview of the differences in resonance due to particle coupling is presented in Figure 18.
Excellent research that has treated this topic can be found in [108,111,112] and is partially shown in Figure 18. It is evident that a significant tuning of the absorption range of plasmonically coupled nanostructures is possible and highly beneficial. Another important consideration are the changes from nonradiative to radiative damping that occur for coupled and uncoupled nanostructures, respectively. Repulsive coupling, in particular, can be postulated to lead to a higher extent of nonradiative decay and a higher heating efficiency. Considerable differences to heating will arise when the nanostructures are supported on solid supports, such as membranes, compared to when they are dispersed in a liquid phase. Plasmonically excited and coupled nanoparticles are vital to achieving controllable and efficient conversion into heat. Specifically, interparticle coupling of the field intensities can shift the far-field and the near-field scattering (the latter becoming important for interparticle distances less than 10 nm for Au nanorods, for example) of the plasmonic absorption considerably [113]. A sophisticated seeded technique of nanoparticle synthesis led to the creation of linked Ag nanospheres with differing numbers of proximal nanoparticles that formed a chain. A clear enhancement in the absorbance to the entire UV-visible range was observed for linked particles as compared to separated particles. A steam generation efficiency of 95% was observed for the best samples, with increases in water temperature of up to 100 °C on irradiation with a Xenon lamp as shown in Figure 19.
When comparing Ag and Au with the less commonly used Pt, Pt has been proven to be more effective than gold nanoshells (AuNSs), a commonly used morphology for PT heating. The loss of the specific morphology of the AuNSs, whereas Pt nanospheres maintained their shape, indicates a better suitability of Pt nanospheres for applications requiring morphological stability at temperatures up to 700 °C [115].
Palladium has a tunable plasmon response in the UV-visible-NIR region with high PT conversion efficiency and PT stability [116], the latter of which originates from its high melting point and the correspondingly higher capacity to withstand high irradiation fluences compared to Ag and Au nanoparticles [117]. This has been confirmed in a study wherein Pd  nano hexagonal sheets preserved their morphology unlike Au and Ag nanoparticles (spheres) when exposed to the same laser fluence. Another significant outcome of this work was the substantial tuning of the LSPR absorption (by more than 250 nm) by changing the edge lengths by a few tens of nanometres [118]. Capping agents, growth directing agents, and surfactants were optimized to achieve homogeneity, size controllability, and absorbance tuning. Although the absorption of these nanostructures were in the NIR region, a tuning towards the UV-visible region was possible through the change in thickness of the sheets [119]. An increase in temperature by 20 °C of an aqueous solution with the dispersed nanosheets and irradiated with a laser confirmed their relevance for PT applications.
The plasmonic properties of copper have been investigated extensively through theoretical approaches, but practical uses for plasmonic applications have been overshadowed by its counterparts (Ag and Au). Copper has been shown to be quite attractive for PT energy conversion, wherein differences in morphology control the light trapping efficiency. Facile one-step laser ablated copper can form hierarchical nanostructures with controlled morphology, depending on the laser scanning rate. Absorbance over a broad range between 200 and 800 nm and with efficiencies over 60% has been reported [21,90].
Similar to the coinage metals copper, silver, and gold, aluminium is a highly promising material for exciting strong surface plasmon resonances in the visible and ultraviolet range, as shown in Figure 15. It is used in plasmonic devices [120] such as resonators [121,122], antennas [123], and biosensors [124]. Aluminium absorbs photons in a wavelength range outside that of conventional plasmonic nanomaterials, which can be exploited for PT applications, such as for absorbing this inaccessible part of the solar spectrum. Al also has a very small scattering cross section, resulting in a much higher PT conversion effect compared to, for example, Au and Ag. Even though Al is plasmonic, the absorption cross section is too low. This shortcoming can be circumvented by using morphologies other than spheres, which can minimize scattering, and through increasing the concentration of the particles [125]. Magnesium nanoparticles have multiple resonances across the UV-vis-IR range [126]. One possible shortcoming for Mg is that they are less lossy than Al and Ga [127], which limits their potential for PT heating. Another shortcoming is the easy oxidation of Mg. Although it does not detrimentally affect the plasmonic properties it reduces considerably the thermal conductivity, leading to a reduced PT performance. The potential of Mg NPs for PT heating, hence, remains unexplored. Of all plasmonic metals considered so far, Na and K have been outliers regarding a study for PT applications [60]. Plasmonically, the very high quality factors of Na and K limit their use for these applications as they are less lossy. Practically, the extremely reactive nature of both metals severly restricts their use.
Tellurium nanoparticles, when engineered appropriately, can have broadband UV-vis plasmonic absorption resonances and are also stable against oxidation, presenting good potential for long term PT applications [35]. They have been proven to lead to higher temperature increases than even perfect absorbers such as TiO x , with an 85% absorbance reported in the region between 300 and 2000 nm. A fundamental analysis of the dielectric function revealed a combination of plasmonic (UV region) and Mie-type multipolar resonances (for particles of sizes from 0.1λ to λ) to be the reason for the broadband optical absorption. The wide size distribution (10 to 300 nm) of laserablated Te nanoparticles allowed for broadband absorption, along with the relatively small bandgap (0.35 eV), which also aids in longer wavelength absorption. Such an optical duality, wherein plasmonic absorption (typically revealed as the most efficient absorption of light) and other resonance modes can be coupled, can prove to be a very effective means of realising the ideal nanomaterial for PT heating [128].
Molybdenum exhibits tunable LSPR between 600 and 1200 nm (for a nanosheet morphology) and has potential in PT studies such as in PT therapy. The lower melting point of Mo, compared to Pd, Au or Ag, is a possible shortcoming regarding PT applications, as laser fluences often lead to thermalization temperatures approaching 1000 K [129]. A comprehensive compilation of the different morphologies of metallic nanostructures explored in PT conversion and the associated data, such as particle dimensions, maximum temperatures/temperature increases achieved, as well as the irradiation means is presented in Table 1. The temperatures in the individual studies have been measured using equipment such as thermal imagers [128,129], thermocouples [130], and IR cameras [131]. Predictive calculations of the final temperatures through calculations [132,133] based on the radiation intensity and the material properties (dielectric constant) have also been performed. Understanding the features of different morphologies, how the materials interact with other substances, and eventually how they might be used in PT applications requires primarily a knowledge of the dimensions of nanoparticles and their effects on SPR and PT performance. Depend on the PT applications, particle morphologies have to be chosen to obtain desired attributes such as absorption range, modality of dispersion, cyclability, and required temperature while also keeping in mind the practical aspects such as cost and ease and repeatability of synthesis. In Tables 1-5 (above and below), we have summarized  the most commonly investigated morphologies for various classes of materials.

Multielemental nanomaterials
The manipulation of morphology leading to significant enhancements in PT properties has been researched through various modalities. One such approach is the fabrication of hollow micro/meso-and nanoporous materials, which have benefits in terms of a lower density (negating sedimentation issues when used for water vaporisation for example) than their solid counterparts, enhanced light absorptivity through internal reflection and increased electron scattering, and the important advantage of having a much higher area of vaporisation (as the pore surfaces also act as vaporisation centres) [169]. Novel nanomaterials, such as mesoporous Au/Ag nanostructures, have been investigated. The bimodal pore size distribution (11 and 3 nm), in addition to enhancing the formation of vapour in the nanoparticle interior and allowing the vapour to diffuse faster to the surface, enhance electron scattering (and hence electron-phonon coupling) by acting as a dense collection of scattering centres, reducing considerably the mean free path of electrons. This has led to a fivefold increase in PT heating (calculated from the ratio of the vapour generation rate to the input energy) compared to a solid Ag/Au core-shell morphology. For example, a Pd@Pt bimetal material for vapour generation is shown in Figure 20 and below in Figure 22. A compilation of the different morphologies and temperature increases for metal-metal and metal-semiconductor nanostructures is presented in Table 2.
The mechanism of self-doping, that is, the manipulation of charge carrier density through creation of vacancies is a crucial  means of realizing plasmonic tuning, not possible with nanomaterials composed of a single element [177]. An important example of a self-doped plasmonic tunable material is copper sulfide, and its PT conversion is shown in Figure 20. Copper sulfide in its pristine form does not exhibit strong plasmonic absorbance due to the overlap with interband transitions [178]. However, removing Cu atoms through various methods, such as annealing, results in a multitude of non-stoichiometric stable chemistries, such as covellite (CuS), anilite (Cu 1.75 S), digenite (Cu 1.8 S), djuleite (Cu 1.96 S), and chalcocite [179] (Cu 2 S), all of which have LSPRs in the NIR due to decreasing amounts of negatively charged copper vacancies (with an associated decrease in electron densities) [180]. PT applications have benefitted from plasmonic Cu 2 S, an example of which are CuS nanoflowers supported on a membrane for solar vapour generation. High PT conversion efficiencies (up to 81.2%) were reported, due to the joint contributions from a broadband absorption as well as dispersion on a membrane on the surface (as compared to a vapour generation system based on a CuS nanofluid). Further efforts on PT heating with Cu 2 S-based materials are highlighted in Table 3.
In a similar multielemental PPT conversion study, porewidened anodic aluminum oxide (AAO) templates as shown in Figure 21 were used to synthesize capillary force-driven selfagglomerated Al 2 O 3 /Au membrane nanomaterials composed of nanowires of this composite. The resulting structures exhibited a multitude of PT heating enhancing factors, namely nanoscale gaps of different sizes, which enable broadband absorption (91% absorption between 400 and 2500 nm due to gaps varying from zero to hundreds of nanometers over a few micrometres), low taper angles (ca. 1°) for sustenance of the plasmon wave throughout the nanowire, thereby avoiding reflection losses, low reflection of incident light (7% between 2.5 and 17 µm), and hydrophilicity, which ensures continuous contact of water during the evaporation process. All of these factors contributed to a high solar thermal conversion efficiency of 57% with an irradiance of 20 kW/m 2 .

Doped metal-oxide plasmonic nanomaterials
Metal oxide nanostructures are suitable for PT applications because of physiochemical stability and increased availability, in addition to prospects of self-doping through the creation of oxygen vacancies. Specifically, certain oxygen-deficient nonstoichiometric metal oxides [186] such as WO 3−x , MoO 3−x and TiO x with a strong absorption across the UV-NIR range with good thermal stability are being researched on increasingly. For example, titanium oxide, with a bandgap energy of 3 eV in the far UV, can be tailored for absorption in the visible region by doping with other transition metals with high efficacy for PT applications [34,[187][188][189]. However, heavily doped or selfdoped TiO 2 in itself is plasmonic, as the self-doping process (due to vacancy creation) leads to excess free electrons due to oxygen desorption. It has been calculated that if free carrier densities exceed 10 20 cm −3 , materials can become plasmonic in the NIR region (the relation between the material and medium dielectric constants notwithstanding). Confirming this postulation, TiO 2 exhibits plasmonic absorption when the ratio between oxygen and Ti atoms reaches approximately 1.67:1. Also, through a broad size distribution of the particles, optical absorption efficiencies up to 90% across the solar spectrum can be reached. The added advantage of TiO 2 being a high-loss material makes it a great candidate for PT heating applications, evidenced by a temperature increase of 50 °C under sunlight of the aqueous solution in which they were dispersed [190].
Another interesting photocatalytic PPT study with an additional layer of SiO 2 is shown in Figure 22. This layer prevented charge hot carrier injection from TiO 2 , ensuring energy dissipation of these carriers as heat and also protecting the AgNPs from oxidation and chemical poisoning. The temperature increase provided by the plasmonic heating improved the hydrogen generation rate from water by fivefold. Since hydrogen generation is a common application in the investigation of the photocatalytic properties of TiO 2 , such combinatorial effects can have benefits such as improving hydrogen yield, minimizing the footprint and cost associated with energy production and desalination, as well as improving the robustness of the materials.
Combining metals with metal oxides can lead to multiple benefits, which has led to extensive studies regarding catalysis, sensing, energy storage, water splitting, and solar cells, to name a few. By extension, the enhanced properties of such composites such as enhanced absorption, broadband absorption, enhanced material stability, and cycling performance have direct relevance for PT applications. One study has investigated the combinatorial effects of AuNPs [192] as a PPT converter, TiO 2 as a photocatalyst as well as an internal light scattering enhancer for desalination, and metal oxide nanostructures for PT applications. The enhancements to the photocatalytic function as well as solar thermal conversion makes such composites attractive for PT applications. A comprehensive summary of different metal-metal oxide composites that have been investigated for PT conversion are presented in the Table 4. 19 Pd@ SiO 2 @ Ag membrane size = 164 nm 6 laser/1 W [191] Pristine Ag nanoparticles had a steam generation efficiency of 60-80% under an irradiation power of 810-930 W/m 2 [203]. There was a considerable enhancement after combination with TiO 2 through the LSPR effect. The efficiency also increased due to an increase in absorbance when the volume fraction of Ag was increased due to an increase in internal scattering within the system [204]. Hence, the plasmon resonance properties of doped metal oxides as well as metal-metal oxide composites may yield new performance benchmarks in PT research [205]. Some of the metal incorporated metal oxide nanoparticles are shown in Figure 23.

Transition-group plasmonic nanomaterials
Transition metal nitrides (TMNs), such as ZrN, HfN, VN, NbN, TaN, CrN, MoN and WN, exhibit plasmonic absorption in the UV-NIR region and have been proposed as alternative materials to the conventional metals [208,209]. Ab initio and molecular modelling calculations predict that group-IV and group-VI TMNs have lower hot carrier lifetimes in comparison to group-V TMNs, indicating a stronger electron-phonon coupling and hence a better potential for PT conversion. Research on PT performance of TiN for example, presented it as an efficient solar absorber and the most efficient one among TiN, ZrN, TaN, HfN, and WN due to the higher absorption cross section across the entire solar spectrum. Chemical stability, high melting point, and visible-to-NIR plasmonic absorption (500 to 1100 nm) of titanium nitride nanoparticles make them attractive, and more so when considering that their scattering efficiencies are comparable to those of Au nanoparticles [210]. Although detailed investigations regarding conversion efficiency, damping properties, and morphology optimisation were not carried out, a threefold temperature increase was observed for aggregated composite TiN particles compared to the control sample without any nanoparticles [211]. An illustration of the PT conversion efficiency comparing the nitrides to Au under the same laser fluence validates their superior properties in this regard, as shown in Figure 24.
Interestingly, the absorption wavelength range of TiN is much higher than that of Au, making TiN highly attractive for PT  Figure 23 panel i was redrawn from [169]. (ii) Sandwich hydrogel with Cu/Carbon cell. Figure 23 panel ii was redrawn from [206]. (iii) TiN on AAO matrix. Figure 23 panel iii was redrawn from [207].  Figure 24 was reprinted with permission from [209], Copyright 2020 American Chemical Society. This content is not subject to CC BY 4.0.
heating applications. This is specifically applicable for PT heating for biological applications, wherein the spectral match between the dipolar resonance of TiN (around 700 nm) and the biological transparency window [212] (650-950 nm) renders it attractive (see Table 5). Detailed investigations of the PT heating of TiN for applications other than for biology are limited. A major limitation of TiN is that it is conventionally synthesized through physical vapour deposition processes such as pulsed laser deposition, reactive sputtering, or molecular beam epitaxy, translating to high equipment and handling cost. This also dictates that they are processed as deposited films, requiring substrates with epitaxial capabilities. Although such substrates are available (Al 2 O 3 and MgO), there is still a restriction on the final configuration of the PT heating system as it must be compatible with films. The availability of solutionprocessed and even doped TiN and other transition group plasmonic nitrides is continually being worked on [212][213][214]. Hence, one can expect new synthesis procedures for highly temperature stable and efficient PT nanomaterials. Combinations of old and new plasmonic materials, each with their own merits in terms of absorption properties, tuning, chemical and thermal stabilities, and broadband absorption, are also being explored, as revealed by studies employing both AuNPs and TiN NPs. The synergistic effect of the resulting composite [197] yielded a better PT heating than either of the pure materials, primarily due to the enhancement of the absorption range of solar radiation. It is evident from the few examples presented here that considerable research effort on plasmonic nanoparticles for PT applications is made, with an unrestricted search on material classes. A detailed compilation of all morphologies and chemistries explored for PT heating [142] is presented in Table 6.
We have graded the performance in PT applications of the materials elaborated in Table 6. Evidently, morphologies other than spheres can offer significant enhancements to PT performance. This can be due to enhancements to the damping of the   resonance (due to scattering with the surface in case of 2D nanostructures/confined nanostructures), enhanced surface area for irradiation, and reduced tendencies for agglomeration. Combined with the effects of plasmonic field coupling on the resonance properties, opportunities abound with respect to the possibilities of employing a combinatorial approach of such morphologies for achieving improved PT performance.

Stability of PPT materials
The definition of nanoparticle stability depends on the desired property that is to be preserved. For catalytic applications, stability pertains to the preservation of the surface area/active crystal facets, whereas for applications such as light harvesting, the active area shares its importance with the necessity of preservation of the optical absorbance properties. Since for PT applications, the latter is more important, stability of nanoparticles for PT applications in terms of morphology and composition is most important, as the optical properties are determined by them, and should be maintained under environments of temperature, pressure, and chemically active atmospheres present in various PT applications [217]. Chemically active atmospheres can, for example, result in a quenching of the plasmonic resonance as a result of strong interactions between adsorbed species and the plasmonic nanoparticle. Anisotropic morphologies such as nanorods [218] or stars [219] can be very susceptible to high temperature, resulting in rounding with a loss of anisotropy and loss of the plasmonic absorbance. Doped plasmonic nanoparticles are particularly vulnerable to temperature and reactive atmosphere effects, due to diffusion, segregation, or reaction of the dopant ions resulting in conversion to a nonplasmonic composition.
More often than not, absorbance tuning of plasmonic nanoparticles is done by effecting a change in the refractive index of the medium, that is, by surrounding the NPs with different materials to form structures such as core-shell [220], hierarchical Figure 25: Colloidal stability of nanoparticles for different particle sizes from 8 to 68 nm of organic pigment samples. Figure 25 was reprinted with permission from [224], Copyright 2017 American Chemical Society. This content is not subject to CC BY 4.0.
[90], or decorated structures [221]. Thus, the stability of the encapsulating material also becomes crucial, as a loss of stability can result in the loss of the tailored absorbance. As an example, highly defect-rich nanoparticles such as twinned/ icosahedral structures [222] are very susceptible to chemical attacks (e.g., by O 2 or salts) and have been reported to be rapidly dissolved, losing all effectiveness for PT heating. Different stability aspects of nanoparticles that affect their optical absorbance will be presented subsequently.

Aggregation stability of nanoparticles in liquidphase PT heating
A common approach to exploit the PT heating effect of plasmonic nanoparticles is to disperse them in PCM hosts, which absorb and store heat during a phase change, such as from solid to liquid or solid to solid. The most common mechanism through which nanoparticles added to liquid PCMs lose their properties is by aggregation ( Figure 25) during storage or during the phase change process. Aggregation stability depends on the collision frequency between the particles and the interaction pair potential between them (a summation of the van der Waal's potential, electrostatic interaction potential, osmotic potential, and the elastic potential), and on whether this potential exceeds a kinetic energy barrier. To obtain this maximum energy barrier, the collision frequency, dictated by the number density of the nanoparticles, and their root mean square velocity ⟨v rms ⟩ is calculated first by [223]: (44) where k B is the Boltzmann constant, T is the temperature, and μ is the reduced mass of two objects.
From the resulting number of collisions per unit time and knowing the desired time over which stability against aggregation is required, the number of collisions during this time can be calculated. Assuming a lower probability of collisions than this over the desired time period will allow for the determination of the threshold barrier for aggregation. Nanoparticle aggregation depends on whether the collisions between the nanoparticles are elastic or inelastic. Elastic collisions result in the colliding nanoparticles remaining separated, while inelastic collisions result in nanoparticle aggregation. Hence, inelastic collisions result when the interaction pair potential exceeds the barrier (V max ) and nanoparticle aggregates are formed [225]. As an example, V max for nanoparticles that need to be non-aggregated for a period of one week has been calculated as [223,226] Thus, the stability of nanoparticles against aggregation can be determined by such calculations, and their handling can be modified accordingly by changing conditions such as pH or by introducing capping ligands/stabilizers [227]. It is also interesting to note that in the case of highly anisotropic nanostructures, such as nanostars, the smallest dimension of the nanostructure in solution determines its agglomeration stability [223].

Chemical stability of nanoparticles in PT heating
It is obvious that the high surface area of nanoparticles leads to a weak chemical stability, which is conventionally viewed as the inherent resistance to irreversible reactions. This resistance is crucial for preventing loss of optical properties. Chemical stability, like agglomeration stability, is influenced by kinetics as well as thermodynamics as shown in Figure 26. Kinetic stability is determined by the ionization energy. As an example, Ag is more easily oxidized than Au due to its lower ionization energy.
Thermodynamic stability of the phase of a nanoparticle is determined by the free energy of formation of its state. Since an exhaustive compilation of the preferred states of nanoparticles is not the focus of this work, only one example of the thermo- dynamic stability of Ag is considered. The oxidation of Ag has a free energy of formation of −11.25 kJ/mol [225]. With a change in nanoparticle size (e.g., for tuning the optical absorbance), the free energy of formation is influenced even at 25 °C [225]: (46) where r is the radius of the nanoparticle in nanometres. Yet another plasmonic metal that is highly promising for ultraviolet PT conversion is aluminium. However, aluminium is easily oxidized with consequent loss of the plasmonic absorbance. This problem is even stronger when aluminium is directly exposed to water because of the significantly enhanced oxidation rate (forming hydrates). Interestingly, the tuning of the nanoparticle dimensions can influence this stability against oxidation. As an example, lithographic aluminium nanodiscs of 100 nm diameter showed a remarkably sustained plasmonic resonance even after a year, compared to larger and smaller discs where the resonance disappeared within days. The reduced grain boundary coarsening [230] of these specifically sized discs due to self-annealing was the reason behind the reduced oxidation rates.
It stands then that nanoparticle chemical stability changes with size/morphology/encapsulation. Nanoparticles can be added with morphologically controllable exterior layers of functional ligands to prevent access of potentially harmful species, such as functionalizing with thiosulfate ions [231], use of conventional surfactants and ligands such as PVP [232], CTAB [233] and citrate [234], and biomolecular ligands such as DNA [235] (due to the strong negative charge on its phosphate backbone) and proteins [236]. Shells composed of organic [237], inorganic [238], or metallic [239] materials are also an option. Another approach to enhance chemical stability is to add dopants such that charge transfer and stabilization through the dopant are achieved. An example is the addition of Ag or Au to stabilize Cu in its pristine [240] form, morphological transformations notwithstanding. Thus, approaches of encapsulation and/or doping can stabilize plasmonic nanoparticles.
The major shortcomings of these approaches that hinder their application in PT heating are the changes to the optical absorbance of the primary nanoparticle as a direct result of the doping or encapsulation processes. These changes include quenching, broadening or narrowing of the absorption, shifts to different wavelengths, as well as changes to other PT-relevant properties such as thermal conductivity, especially in the case of the encapsulation. Changes to thermal conductivity due to encapsulation can be a bigger concern than doping, as loss of thermal conductivity due to phonon scattering at the interface of the encapsulating layer will be present for these materials, but not necessarily after doping. A careful selection of the appropriate methodology for enhancing the chemical stability is therefore essential.

Thermal stability of nanoparticles for PT applications
As discussed in the previous sections, the nanoparticle morphology and the environment control the optical absorbance. Any changes to these factors hence can be detrimental to the suit-  Figure 27 panel ii was reprinted with permission from [244], Copyright 2020 American Chemical Society. This content is not subject to CC BY 4.0.
ability of these particles for the concerned purpose. The change in surface energy is a driving force for shape changes in nanoparticles. Despite preservation of the desired optical properties due to controlled synthesis, post-synthesis shape changes can happen through agglomeration (during storage or prolonged use), chemical reaction, and (often more frequently) thermal factors such as increasing temperature. Indeed, the latter mechanism has been developed as a nanomaterial fabrication technique on its own, based on the PT effect and termed as light-induced rapid annealing (LIRA [142,241]). A well-established mechanism of particle size increase due to temperature is the phenomenon of Ostwald ripening, wherein atoms on the surfaces of smaller particles dissolve and deposit on the more stable surfaces of larger particles. This mechanism has not been conclusively proven for nanoparticles dispersed in a liquid phase, although particle growth due to agglomeration presents a significant challenge. For nanoparticles dispersed in a solid matrix (such as the case of solid PCMs), this process can be inevitable without the presence of stabilizing agents/encapsulating layers. Examples of this phenomenon include Pd nanoparticles ripening at ambient temperature [242] as well as at high temperatures of 750 K [243]. Heat trapping layers can be added to achieve plasmonically enhanced PT efficiency and to increase the thermal stability of the plasmonic materials, while also reducing oxidation. For example, Al 2 O 3 has been added as a trapping layer on Ag to increase the stability to over three months, while Ag alone underwent oxidation within 20 days [198]. The thermal stability of Au/Ag nanoparticles with different morphologies is shown in Figure 27 and Figure 28.
Temperature-driven interdiffusion influences significantly the thermal stability of a nanomaterial, especially when present as alloys or as doped materials. The Kirkendall effect [247] dictates that the interface between two materials having different diffusion rates moves as a function of the temperature as well as the nature of the environment. The result is the formation of new chemistries as well as morphologies. The Kirkendall effect can be envisioned to be a positive or a negative mechanism, wherein the positive aspect is the possibility to achieve novel morphologies such as nanoboxes, fullerene-like nanoparticles [248], and porous structures [249] through facile processes. The negative aspect then is the loss of PT properties due to changes in chemistry and/or morphology even at room temperature. Hence, the knowledge of the diffusion rates of materials is important when employing multiple compositions/ phases during the design of PT plasmonic nanomaterials. An evident, but often infeasible, solution to avoid interdiffusion is the storage and handling of the prepared nanomaterials at low temperatures. (ii) (a-c) To prevent this, AuNPs can be incorporated into support materials such as yttrium-stabilized zirconia, titania, or resorcinol formaldehyde (RF). Figure 28 panel ii a was adapted with permission from [114], Copyright 2019 American Chemical Society. This content is not subject to CC BY 4.0. Figure 28 panel ii b was adapted with permission from [246], Copyright 2017 American Chemical Society. This content is not subject to CC BY 4.0. Figure 28 panel ii c was adapted with permission of SPIE from [218] ("Investigation of the optical and sensing characteristics of nanoparticle arrays for high temperature applications", by G. Dharmalingam  Thermal stability of a nanomaterial in PT heating depends on the melting point as well. The diffusion of atoms becomes unrestricted on reaching this temperature at any location within the nanoparticle. Melting is often explained as liquefaction when the amplitude of atomic vibrations exceeds a fraction of the interatomic spacing. The influence of temperature on this amplitude and hence on melting was established as early as 1910 by the relation [250] where m is the atomic mass, ν E is the Einstein frequency, δx is the root-mean-square thermal average amplitude of vibration, k B is the Boltzmann constant, and T is the absolute temperature. For nanomaterials, the melting temperature is closely associated with the particle size and given by the relation [251] where T m and T mbulk are, respectively, the melting temperature of the nanoparticles and bulk material, D is the diameter of the nanoparticles, and C is a parameter that depends on the model applied for the melting process, that is, the homogeneous melting hypothesis (HMH), the liquid skin melting (LSM) model, or the liquid nucleation and growth (LNG) model [252]. This Gibbs-Thompson relation has been used successfully in multiple studies for estimating the depression in melting point with the size of nanomaterials. It is to be noted that although this model is considered sufficiently accurate by multiple studies, variations of modelling the melting temperature exist, depending on how the melting point is arrived at, such as by modelling the Lennard-Jones potential [251], or whether the microscopic variations in thermal properties and aspects such as shape effects [91,253] are included in the calculations. An example of alternative studies that have successfully modelled nanoparticle melting is the study of the effect of laser irradiation on Au nanorods, wherein the model could predict the temperatures that cause reshaping of the nanorods (due to plasmonic heating) [143]. Temperatures above 1000 K at the irradiation location reshaped the nanorods, while the lower temperatures at the boundaries of the beam spot (795 K) did not affect them [14]. In a representative study, the significantly higher melting point of palladium nanoparticles led to a better performance and maintenance of morphology (and hence the LSPR property) compared to Au and Ag nanoparticles under radiation with NIR light, as detailed in prior sections [119]. The water evaporation rate of the confined linked Ag-Au nanoparticles remained stable after ten cycles while the morphology also remained stable, demonstrating the potential for encapsulation in enhancing nanoparticle stability for PT applications [114].
Platinum nanoparticles normally absorb in the UV-visible range with extensions into the NIR region [115]. Observations such as the melting of the adjacent areas of a NIR-irradiated Pt nanoparticle at temperatures greater than 700 °C and the stability of the nanoparticle itself leads to the postulation that they can be viable candidates for PT heating applications. While the process of laser irradiance for studying plasmonic heating is the rational approach to obtain precise information on the absorbance properties, it must be kept in mind that an extremely localized heating of the nanoparticles or the nanoparticle dispersion occurs in these studies. This observation, although accurate for applications such as photodynamic therapy, does not reflect accurately the processes such as water vaporisation through solar power irradiation, wherein the irradiation area is considerable higher than a micro-/millimetre-sized laser beam. This difference in the heating process is important as it can lead to very different outcomes. Nanoparticles to be used for solar heating studied by laser irradiation will not necessarily melt and undergo morphology changes when eventually deployed in practice. Hence, it is crucial that nanomaterials to be employed for PT conversion applications are characterized by accurately reproducing the field deployment guidelines.
The discussion of the stability of nanomaterials for PT energy conversion thus takes into account the predominant reasons affecting the same, namely temperature, chemical reactions, and aggregation. Regardless of the phase of the PT conversion system (liquid/solid/gas), these three phenomena govern the suitability as well as the life cycle of the nanomaterial. Proper guidelines that can be derived from the application need to be considered for the selection and design of PPT nanostructures for optimal performance.

Conclusion and Outlook
Plasmonic nanoparticles are definitely a topic in PT research, be it in medicine for therapy or regarding energy for storage, heating, or chemical syntheses. We have presented a review on established as well as emerging plasmonic materials that are of importance for a multitude of applications. Among the principal requirements for a material for PT research is the conversion efficiency. It is clear that among the materials presented, the conversion efficiency of pristine materials is markedly lower than that of composited materials, for reasons such as non-ideal dielectric properties (for heating), performance degradation over time, limited possibilities in terms of the absorption wavelength range, and stability against aggregation, reaction, or temperature. Thus, materials composed of multiple elements and/or morphologies are needed to advance the practicality of PPT materials. Novel morphologies such as linked nanoparticles, collapsed nanowires, and concentric rings composed of multiple materials or even the simple polydispersed nanoparticles can lead to vastly better PT properties in terms of enhancing the absorption range and improving the conversion efficiency. Interesting properties such as optical duality observed in, for example, Te can be exploited to achieve broadband absorption of light through the manifestation of both plasmonic and Mie resonances. We have reported on a narrow hybrid resonance mode, which arises due to coupling of the plasmonic and a second diffuse mode on an underlying metal oxide layer. Thus, hybrid (or) multiple-resonance materials can be promising for PT research. A minor shortcoming of such combinations is the inadequacy of simulation methods to comprehensively predict the optical properties. Often chemical reactions that might occur when bringing such materials together are not known and cannot be included in the simulations. Nevertheless, simulations can help in understanding discrepancies between expected and observed material behaviour. It is the authors' opinion that a major hindrance to broadening the spectrum through combining materials for achieving optimal plasmonics is the stability or lack thereof when these materials are brought together. Deep understandings of the aspects which lead to instabilities, that is, aggregation, irreversible reaction, and thermal degradation can help in this regard. This is also valid for the unexpected interplay between grain size and the enhancement in PT property preservation. With continual pursuits on improving the presented aspects, plasmonic materials can be expected to drive PT research towards a better future.